Single point potentials with total resonant tunneling 
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Two rectangular models described by the one-dimensional Schrodinger equation with sharply 
localized potentials are suggested. The potentials have a multi-layer thin structure being composed 
from adjacent barriers and wells. Their peculiar tunneling properties are studied in considerable 
detail. Particularly, in the zero-range limit when the potentials are squeezed to a single point, 
sharp peaks with total transmission are observed at certain (positive and negative) quantized values 
of the potential strength constant forming infinite discrete sets. Beyond these sets, the barrier- 
well structures behave as a perfectly reflecting wall. The transcendental equations with respect to 
potential strengths, the solutions of which determine transmission (resonance) sets, are derived. In 
this regard, both the models are exactly solvable. The energy dependence of an incident particle 
is shown to reveal a resonance behavior, being completely different from that observed in a typical 
double barrier structure. 
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I. INTRODUCTION 

The mechanism of resonant or non-resonant tunnel- 
ing through two or more successive potential barriers 
plays an important role in studying different phenom- 
ena in quantum physics and optics including various as- 
pects of engineering at micron and nanoscale levels 
Oriented towards the investigation of the resonant tun- 
nel effect, beyond the standard case of double barrier 
(DB) penetration, it is of interest to search for the mod- 
els where other mechanisms of resonant transmission can 
take place. Thus, a one-dimensional model of resonances 
described by the Schrodinger equation with Dirac's delta 
function barrier and an abrupt mass jump has been sug- 
gested and studied in Ref. Q. It is worth mentioning 
that the motion of a particle with a position-dependent 
mass under different conditions (for instance, including 
an oscillating field) with further applications to nanos- 
tructures was a main object of studies in recent publica- 
tions H. Another example of recent studies relevant to 
resonant tunneling are one-dimensional sharply localized 
potentials consisting of barriers with adjacent wells (3~ 
6]. Here the approach is based on the search for those 
finite-range potentials which in the zero-range limit ex- 
hibit peculiar features like a non-zero transmission under 
certain conditions. 

In general, the models with zero-range potentials that 
describe point or contact interactions (for details and ref- 
erences see books @, HI), are widely used in various ap- 
plications to atomic physics condensed matter sys- 
tems liql. acoustic and electromagnetic waveguide net- 
works [11]. Currently, there is a growing interest in the 
quantum dynamics of particles in structures of low di- 
mension which can be described by graphs (for details 
see survey fl2j|) using different co uplings and their ap- 
proximations at the graph vertices |13l4l5l | . Of particular 
importance are nanoscale quantum devices such as thin 
quantum waveguides [l(| and potential-controllable spec- 
tral filters [T3. 



The point interaction models admit exact closed ana- 
lytical solutions, which are rare but quite useful. They 
provide relative simple situations, where an appropriate 
regularization procedure can be chosen to be in rele- 
vance with real structures being important for technolog- 
ical applications. Thus, remarkable structural analogies 
between the Dirac delta function and Coulomb poten- 
tials in one dimension have been discussed (18J. In the 
case of moving this delta potential with a constant ve- 
locity, the exact propagator has been evaluated by sum- 
ming over the spectrum of eigenstates and by calculating 
the path integral directly [19| . In a number of papers 
poj, the relativistic case of a delta function potential 
in the one-dimensional Dirac equation has been treated. 
On the other hand, for a general class of point inter- 
actions, Green's functions have been computed by us- 
ing the reflection and transmission coefficients 21]. The 
Schrodinger equation with two- and three-dimensional 
delta function potentials has been the object of consider- 
able study [22]. In particular, exactly solvable models has 
been elaborated in two and three dimensions by renor- 
malizing the kinetic energy operator [23| . It is also worth 
to mentioning the recent studies of time aspects of point 
interactions [24j . 

In this paper, we would like to address the problem of 
resonant tunneling through single point potentials. The 
regularized (finite-range and realistic) structure of these 
potentials has the form of barriers with adjacent one or 
more wells. The presence, at least, one well will be shown 
below to play a crucial role in the mechanism of tunnel- 
ing. The key point of our approach is the understand- 
ing of this mechanism which appears to differ from that 
well known in DB structures. A small region of large re- 
pulsive potential (barrier), which is followed by a small 
region of large attractive potential (well), can be served 
as an interesting example in one dimension. Under an 
appropriate parametrization of parameters, this barrier- 
well (BW) potential in the zero-range limit becomes the 
derivative of Dirac's delta function referred to as a point 
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dipole studied in numerous publications [25|, |26(. An- 
other example, which can be treated analytically, is the 
second derivative of a delta function but with a renormal- 
ized strength constant @, which essentially reduces the 
singularity at the origin of the point interaction. Below, 
it will also be demonstrated that the zero-range idealiza- 
tion is an useful prompting for investigation of tunneling 
through finite-range structures. 

Thus, we consider the systems that present two suc- 
cessive barriers with a well inserted between them or a 
barrier with adjacent wells located at both its sides. More 
precisely, these systems are supposed to be described by 
the following two potentials (see Fig. [TJ: 



A+(x) = A' e (x + r) + A' e (x-l) 
A-(x) = A' e (x + r) + A' s (r-x), 



(1) 



given through the rectangular function A' e (x) of the BW 
form defined by the equations 



h for - I < x < 0, 
A' s (x) = { -d for < x < r, 
otherwise. 



(2) 



Here, h and I are height and width of the barrier, and 
d and r are depth and width of the well, respectively. 
These parameters can be expressed through a squeezing 
parameter e as follows 



h = 
d = 



2^ 



cx{c\ + c 2 ) 

2a- 
ca(ci + c 2 ) 



I = C±E, 



r = c 2 e, 



(3) 



with arbitrary positive constants c\ and c 2 . The con- 
stants a^- are defined as follows: (a + , er~) = (1, 1) for 
a = 0, (er + , <j~) = (1, a) for a > 0, and (<r + , er~) = 
(<r, 1) for a < 0, where a is an arbitrary positive pa- 
rameter. This definition means that the parameter a is 
assigned only to the wells. Parametrization ([3|) has been 
chosen in such a way that in the particular case a = 1, 
one obtains the weak limit: A' s (x) — > S'(x) = d5{x)/dx 
in the sense of distribution theory as e — > 0. In this 
case, Ih = rd, i.e., the condition j_ A' E (x)dx = is 
fulfilled. The parameters (or a) have been incor- 
porated to investigate a more general class of point in- 
teractions than the <5'-potential. They control the devi- 
ation of the standard §' (x) distribution from the zero- 
range limit of the whole cr-family of point interactions 
which cannot be defined as distributions. Furthermore, 
the parameter a appears to be useful in studying the 
situation, when the wells continuously disappear result- 
ing in a DB structure as a — > 0. In the particular case 
cr = l, within parametrization ([3]), one can prove that 
A±(x) -> S'(x+0)±S'(x-0) in the sense of distributions 
defined on the space of appropriately chosen test func- 
tions. Therefore, if a = 1, the potentials Af(x) may be 
called double (5'-like structures or one-dimensional double 
point dipoles (or double plane capacitors). 
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FIG. 1: Rectangular models: Af(x) (top) and A £ (x) (bot- 
tom) defined according to Eqs. ([1} and ([2]), and parametrized 
by Eqs. © with cr + = er - = 1. 



Thus, we consider the one-dimensional Schrddinger 
equation for stationary states, which in the limit that 
neglects the interaction between electrons, reads 



%l>" (x) + aAf (x)-4>(x) =E*P(x), 



(4) 



where ip(x) is the wave function for a particle of mass m 
and energy E (we use units in which h 2 /2m = 1), and a 
is an interaction strength constant. 

Recently, using various regularization schemes 
for the distribution S'(x), the potential aS'(x) inserted in 
Eq. (U]) instead of aAf(x), has been shown to be partially 
transparent at a certain set of quantized values for a, but 
acting as a perfect reflecting wall for the values different 
from this transmission set. This resonant transmission 
property has rigorously been proved in Ref. 28] by us- 
ing a general class of regularizing sequences. However, 
the partial transparency decreases very rapidly with the 
growth of a. Consequently, from a practical point of 
view, it would be of interest to search the single point 
interaction systems that admit a total resonant transmis- 
sion. A composition, at least, of two adjacent (^'-function 
potentials is expected to permit a total resonant behav- 
ior. That is why the two double <5'-like structures de- 
scribed by Eqs. (JTJ) and plotted in Fig. [T]have been cho- 
sen for this purpose. In the zero-range limit, for one of 
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these structures the two point dipoles are arranged for- 
ward (top of Fig. [TJ) , whereas for the other one they are 
arranged oppositely (bottom of Fig. Q} . 



II. FINITE-RANGE SOLUTIONS FOR THE 
A^(as) MODELS 

The positive energy solution of the Schrodinger equa- 
tion Q with potentials (fTJ) can be written through the 
transfer matrices A^, connecting the boundary condi- 
tions for the wave function ip(x) and its derivative ip'( x ) 
at x = x\ = — I — r and x — X2 = I + r: 



ip(x 2 ) \ _ 



= A= 



ip{xi) 
I/fa) 



A 



= ( A | 

A 21 ^22 



.(5) 



This approach appears to be very useful because the re- 
flection and transmission coefficients can be computed 
immediately in terms of the transfer matrix. On the 
other hand, the zero-range limit of transfer matrices di- 
rectly leads to the corresponding connection matrix, on 
the basis of which all possible (four-parameter) point in- 
teractions in one dimension can be classified j29|. Simi- 
larly, the scattering data can easily be calculated thr oug h 
the connection matrix (see, for instance, Refs. [5l. l2ll . l26l ]). 

Thus, by matching the values and derivatives of the 
wave function i^(x) at each potential discontinuity, after 
straightforward calculations we arrive at 



cos(2pl) cos 2 (qr) 



5 V <7 + P ' sm,2 ^ /!su ' ,2 ' / '' ) 



sin 2 (pi) — cos 2 (pi) 



sin 2 (qr), 



\ + 

A 12 



- sin(2pl) cos 2 (qr) + - cos 2 (pi) sin(2gr) 



■ sin(pZ) cos(qr) 



(5- 






p) 



1 



cos(pl) sin(gr) 



svet(pl) sin(gr), 



A 21 = — psmi2pl) cos 2 (qr) — qcos 2 (pl) sin(2qr) 
— I — ] \psm(pl) cos(qr) 

q pj 



A 22 



qcos(pl) sin(gr)] s'm(pl) sin(gr), 
cos(2p^) cos 2 (gr) 

2 (f + 3 p) sin<y2pl ^ sin ( 2, ? r ) 



sin 2 (pi) - cos 2 (pi) 



sin 2 (qr) 



(6) 



and 



A u = A 22 = cos(2pl) cos(2qr) 

_ I (? + £)sin(2pl)sin(2gr), 
2 \q P) 



Aio — 



A 21 — 



- svt\{2pl) cos(2qr) 
P 



-cos 2 (pi)- - sin 2 (pi) 
q p 



—p sin(2pl) cos(2qr) 
2/„n 1 „„„a 



P 



— sin 2 (pi) — - cos^(pZ) 
q p 



sin(2qr), 



sm(2qr), (7) 



with 



P 



\[E - ah , q = y/E + ad , E = k 2 



(8) 



Here, p and q can be either real or imaginary depending 
on the sign and values of a and the wave number k. 

According to parametrization (J3J , the quantities p and 
q are divergent as e — ¥ 0, while the products pi and qr 
must be constant in this limit. Under the latter assump- 
tion, it follows from the form of Eqs. ^ and (J7J that 
the matrix elements A^ and \ 22 are constant, X^ 2 0, 
and A^ is divergent as e — > 0. However, one can expect 
that under some appropriate conditions, the cancellation 
of divergences may happen. Indeed, the terms A^ can 
be factorized as 



A 



[(!. 




\q 


pj 



sin(pZ) sin(qr) — 2 cos(pl) cos(qr) 



x \psm(pl) cos(qr) + qcos(pl) sin(gr)] , (9) 



— sinfpZ) sin(gr) — cos(pl) cos(qr) 

q 



A 21 — 2 

x [p sm(pl) cos(qr) + q cos(pl) sin(qr) 



(10) 



Using next this representation and imposing the con- 
straint A^ = 0, we find the conditions under which the 
divergences cancel in the e — > limit. As a result, for 
the potential aA+(x) one can write the following two 
equations: 



- + - I shxOpl) sin(gr), (11) 

q p; 



2 cos(pl) cos(qr) 
p sm(pl) cos(qr) = —q cos(pl) sin(qr). 



(12) 



Similarly, for the potential aA £ (x), the constraint A 21 = 
also leads to the two equations: the same Eq. (fT2")) and 



p sin(pZ) sin(gr) = q cos(pl) cos(qr). 



(13) 



Next, imposing the condition A^ = and using 
Eqs. (rri| - (fT3"|) . one can transform the rest of the matrix 
elements in Eqs. © and ([7]), yielding the representations 



At 



A 2 + = 



a; 



cos (pi) 
cos 2 (qr) 



-A" 



cos (qr) 
p sin(2p() 



cos (pi) 



-l 
o 



q sin(2gr) 
cos 2 (qr) 
cos 2 (pi) 

2p~ 1 (p 2 -q 2 ) tan(pl) 
' p 2 tan 2 (pl)-\-q 2 " 

- 1 



(14) 
,(15) 

• (16) 
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where I is the unit matrix. Here the representations 
and A 2 for the aAf(x) potential have been obtained 
under conditions (fTTj) and (fT2]l . respectively. Similarly, 
the representations and A^ for the aA~ (x) potential 
correspond to Eqs. (fl~3|) and (fT2)) . respectively. 



III. 



TRANSMISSION SETS IN THE 
ZERO-RANGE LIMIT 



In this section, we analyze the e — > limit of Eqs. (fTTj)- 
(fl~3| as well as A^ 2 -matrices (fT5]) and (fT6|) . As a result, 
using the asymptotics 



p — >• 1 , 



ci(ci + c 2 ) 



£ 1 and g 



2aer 



C2(Ci + c 2 ) 



(17) 

for the potential aAf(x) with any a g (—00, 00), one 
finds that the £->0 limit of Eq. (fTTj) yields the equation 



ba~ 

(7+ 




tanh 



2ao-+ /2ao- 
— tan-i 



1 



1 



(18) 

where 6 — c\jci. Similarly, for the potential aA e (x), we 
get the e -> limit of Eq. ffl3]) : 



tanh 



2a<r J - 



tan 



2acr" 



Equations (1181) and (|T9|) admit an infinite but countable 
set of solutions (roots), which we denote by (n — 
0, ±1, . . .) for the potentials aAf(x), respectively. In its 
turn, Eq. (fT2j) also determines a countable set of quan- 
tized values for a, being the same for both the potentials 
aA^(x). The corresponding limiting equation becomes 



tanh 



2ckj+ 
1 + 6- 



ba / 2aa 

— — tanW - 

cr+ V 1 + 6 



(20) 



We denote this set, which in the particular case a + = 
<j~ = 1 coincides with the corresponding resonance set 
found in Ref. [27] for the <5'-potential, by a' n (n — 
0, ±1, . . .). However, the zero-range limits of the A 2 - 



matrices described by Eqs. (jT5j) and (1161) are different. 
Since the off-diagonal matrix elements \f 2 in these equa- 
tions disappear as e — > 0, the connection matrices A^ 2 
take the following final form: 



lim A? 



-I, limA+ = A', limAj=7, (21) 



£^0 1 E-»0 £^0 

where the A'-matrix is given by 

, 



A' = 







^ cosh^2<a+/(l + 6- 1 ) 1 



cos^/2o4o-/(l + 6) 



(22) 

In the particular case ct + = <t =1, this connec- 
tion matrix corresponds to the ^'-potential computed in 
Refs. @,H3,[2|. 



Note that the existence of Eqs. (|T8|) - ([20)) . which deter- 
mine the conditions for a non-zero transmission, imme- 
diately follows from the form of the potentials a Af (x) . 
Indeed, due to the scaling given by Eqs. ([3]), in the zero- 
energy resonant case E — > 0, Eq. ([5]) becomes homo- 
geneous and scale invariant with respect to squeezing 
the potentials. Therefore for the existence of non-trivial 
solutions, it requires a constraint on the potential pa- 
rameters a, c\, C2, and ct* fixing admissible sets for a. 
In our case, the constraint can be found explicitly by 
matching the rescaled wave functions and its derivatives 
at x/e = £ = 0, ci, ±C2, ±(ci + C2). As expected, within 
this straightforward procedure we arrive at Eqs. I|18p- 
420]). 

Thus, as follows from Eqs. (|2Tj) and (|2"2")h the point 
potentials determined by the zero-range limit of the 
Schrodinger equation (UJ) describe the total resonant 
transparency on the sets E ± = { a n}^—-oo f° r both the 
potentials aA^(x), respectively, and on the set E' = 
{o4}$?L-oo f° r the potential aA~(x), whereas only the 
partial transparency occurs for the potential a Af (x) on 
the set E'. Beyond these transmission (resonance) sets, 
A 21 —> 00 as e — > and therefore Eqs. (|lip - (fT5|) are not 
valid anymore. For these continuous values of a (almost 
everywhere on the line — 00 < a < 00, except for the 
resonance sets E ± and E'), both these point potentials 
(referred to as "separated" while using either the local 
[29 ] or the global [lj] description of the boundary condi- 
tions at the singularity point) are opaque. 

Note that A'-matrix (|22p connects the two-sided 
boundary conditions at x = ±0 as 

V->(+0) = 0„V(-O) and V'(-0) =M>'(+0), (23) 

which mean that the wave function and its derivative are 
discontinuous at x = 0. As described by these equations, 
the discontinuity is expressed in the special way through 
the function 9 n defined only on the discrete set E'. In 
more general terms, for the one-dimensional Schrodinger 
equation, the approximation of point interactions with a 
discontinuity of wave functions at the singularity point 
has been analyzed earlier in Ref. (30j . where instead of 
rectangles, the sum of three delta functions: 

£(x; a, c, e) = c(e)5(x + e) + a(e)S{x) + c(e)5(x - e) (24) 

has been used. Within this approach, varying the behav- 
ior of a(e) and c{e) as functions of e in a proper way, a 
whole class of point interactions, which realize a discon- 
tinuity of wave functions as the distance e between the 
delta functions tends to zero, has been analyzed. 



IV. THE SCATTERING APPROACH: 
NUMERICAL CONSIDERATIONS 

For better understanding the phenomenon of tunnel- 
ing across the aA^ (x) potentials, it is reasonable to carry 
out direct numerical computations of the scattering am- 
plitudes when the parameter e is finite. The reflection 
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and transmission coefficients can easiiy be expressed in 
terms of the elements of the connection matrix for a given 
point interaction [U [2(| . Similarly, for finite-range in- 
teractions, these coefficients can be computed through 
the transfer matrix; in our case, by employing the A^- 
matrices described by Eqs. <(6|)-((8]). Thus, omitting for a 
while the superscripts " ± " , we consider the scattering 
both from the left (labelled by the subscript "/" ) and from 
the right (labelled by the subscript "r" ) of the structures 
described by Fig. [TJ As usual, we adopt the following 
definition for the reflection and transmission coefficients: 



ip(x) 



Tie' 



Ri e 



-ikx 



for 
for 



— OO < X < X\, 
X2 < X < OO 



from the left and 

-ikx 



ip(x) 



c 
T r e 



+ R r c lkx for X2 < x < oo, 
lkx for — oo < x < xi 



(25) 



(26) 



from the right, where x\ — —l—r and x-i — l+r. Inserting 
Eqs. (f25|) and (|26| into A- matrix equation ([5]), we find 
the reflection (Ri and R r ) and transmission (T/ and T r ) 
coefficients: 



Ri = D- 1 [A 22 -A 11 -i(fcA 12 + fc- 1 A 2 i)]e : 
D- 1 [A n - A 22 - i (AAia + fc^Aai)] 
T r = 2D- 1 c ik{xi - X2) , 



R r 
T, 



2ikx± 
-2\kx2 



(27) 



where D — An + A 22 — i (A;Ai2 — fc _1 A 2 i). Using now the 
relation AnA 22 — Ai 2 A 2 i = 1, the probabilities of reflec- 
tion (reflectivity TV) and transmission (transmissivity T) 
can be written in the form 



U = \R,V = \Rr\ = 



t = \TA 2 = \T r \ 2 = 



u 2 + v 2 

+ U 2 + I 

4 



where 



An — A 2 



and 



feAi 



(28) 



(29) 



The equation for the transmissivity T, where the A^- 
matrices are given by Eqs. ©-(HI), will be used below for 
numerical computations by varying both the parameters 
a and k. 

Now we focus on the direct numerical calculation of 
the transmissivity T — T ± according to Eqs. (|28|) and 
(|2"9")l . The quantities u — u ± and v ~ v are computed 
by using Eqs. <j3j) and ©-(IH]). For better visualization, 
the logarithmic scale has been used. The case with a = 1 
provides the double (5'-like structure, while the case with 
(7 = is relevant to the absence of wells. 

First, at fixed values of b and k, we calculate the trans- 
missivity T as a function of a for some values of e > 0. 
Thus, T + plotted in Fig. Has a function of a for two 
fixed values of e clearly illustrates the five peaks with 
total transmission (T + = 1) and the three peaks with 
partial transmission (T + < 1). These peaks correspond 




FIG. 2: (Color online) Transmission probability T + plotted 
as a function of the strength a at fixed wave number (k — 1) 
and e = 0.1 (curve 1, black) and e = 0.2 (curve 2, red). Five 
peaks with total transmission (T + = 1) at al 2 — —18.26 
(inset, bottom), «1 1 = -2.70, = 0, a+ = 2.28, a\ = 
35.09 (insert, top) and three peaks with partial transmission 
(T + < 1) at q'_ 2 = -35.82, a!_ Y = -11.66, a[ = 26.87 are 
clearly visible. Here 6 = 3 and a = 1. 



to the solutions of Eqs. ((T8)) and (|20|) approaching them 
in the e — > limit. In the natural (non-logarithmic) 
scale, these peaks appear to be very narrow even if e is 
not small. The location of the peaks exactly corresponds 
to the sets E + and In the space between the peaks, 
the transmission is almost zero. 

It is unusual that the transmission equations (fT51) - (|2Tj)) 
do not contain any dependence on the energy of an in- 
cident particle E (or the wave number k). Therefore we 
have also studied the behavior of the transmission as a 
function of k in a realistic case when e is not very small. 
To this end, we have plotted the dependence of the trans- 
mission on both the quantities a and k. In the form of 
contour plots, this dependence is demonstrated by Figs. [3] 
and @] for both the potentials aA^(i). These figures, 
where the ridges show the highest transparency, illustrate 
the similarity of the tunneling behavior. At small k, the 
ridges are located exactly on the sets E + and £' (Fig. [3]), 
and E~ and £' (Fig. 2]). As further demonstrated by 
Figs. [3] and 21 for low particle energies, any dependence 
on k is absent, but with increasing k, the ridges shift 
to smaller a compared with the solutions to Eqs. (fi"8|)- 
(l2Cfl) . The central region corresponds to overbarrier trans- 
mission which is obviously transparent. The subbarrier 
tunneling is found in the region bounded above by the 
curves: 



k < 



2n 



ci(ci+c 2 ) 



2|q| 



for a > 0, 
for a < 0. 



(30) 



e 2 (e 1 +c 2 ) 

In order to illustrate the profile of the ridges described 
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FIG. 3: (Color online) Transmission probability T + plotted 
as a function of the strength a and the wave number k with 
e = 0.2. Here b — 3 and a = 1. The ridges at k —¥ exactly 
correspond to the peaks shown in Fig. [2] 
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FIG. 4: (Color online) Transmission probability T~ plotted 
as a function of the strength a and the wave number k with 
e = 0.2, 6 = 3, and a = 1. The ridges at — > correspond 
to total transmission and they are given by the solutions to 
Eqs. (HU and (2D]): a'_ 2 = -35.82, aT 2 = -11.74, aLj = 
-11.66, aZi = -1.01, a = 0, a7 = 8.77 and c*i = 26.87. 



by Figs. El and SI we have plotted the dependence of the 
amplitude T + on the wave number k at some fixed val- 
ues of a being close to the solutions of Eqs. ([T5]) and 
([20]) . We have restricted ourselves to the case with to- 
tal transmission and found that the corresponding peaks 
are very narrow in the natural scale. However, in the 
logarithmic scale, the peaks essentially spread and the 
profiles become much more convenient for visualization 
as demonstrated by Fig. [5] The peaks in this figure cor- 
respond to the ridges depicted in Fig. |3] Note that these 
peaks correspond to the quantized values of strengths: 
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FIG. 5: (Color online) Three sections of surface T + (a, k) at 
fixed strengths: a = 25 (curve 1, black), a = —16 (curve 2, 
red) and a = —34 (curve 3, blue). The parameters chosen for 
calculations are the same as for Fig. [3] 




FIG. 6: (Color online) Transmission probability 7 plotted 
as a function of the strength a and the wave number k with 
a = 0. The rest of parameters chosen for calculations are the 
same as for Fig. [4] 



a = a^ 3 , a^ 2 , a J at which the total transmission takes 
place. 

It is of interest to clarify the origin of the resonances 
described above for the potentials Ajf (x) with a = 1. To 
this end, we have compared these resonances with the 
well known phenomenon of the tunneling through DB 
systems. The DB structure can immediately be obtained 
from both the potentials aA^(ir) if we remove the wells, 
i.e., put (j = in Eqs. ([B])-©. For instance, in the case 
of the potential aA~ (x) where = AJ 2l the require- 
ment v~ = k\^ 2 + k~ 1 \2i = results in the well-known 
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resonance condition 

(t + -) tanfpZ) = 2cot(2fer). (31) 

V fc p/ 

On the other hand, as follows from the form of 
Eqs. (IT51) - (|2T)1) . all the solutions to these equations run 
along the a-axis to infinity and they obviously cannot be 
obtained from Eq. (|3"Tj) . Consequently, one can conclude 
that the presence of wells adjacent to the barriers makes 
these combined structures transparent at some quantized 
values of the strength a. In this regard, we can refer this 
type of transmission to as a BW resonant tunneling. One 
can suppose that this tunneling differs in its nature from 
the double barrier resonant transmission. This differ- 
ence can be better understood if we compare the contour 
plots given by Figs. U (cr — 1) and [5] (<J = 0), where the 
rest of the parameters are the same. As shown in these 
figures, the ridges of resonant tunneling are "perpendic- 
ularly" arranged with respect to each other. Under the 
well presence, the resonances occur nearby the quantized 
values of the parameter a being very sensitive while vary- 
ing this parameter. In the case of the well absence, the 
a-dependence of the resonances disappears at all. 

V. CONCLUSIONS AND DISCUSSION 

Although rectangular systems are mainly straightfor- 
ward exercises in textbooks on quantum mechanics, nev- 
ertheless, their various parametrizations leading to one- 
dimensional point or contact interactions provide rich 
objects with non-trivial properties. In this paper, we 
have introduced a new type of single point interactions 
with a fully resonant transparency. Their finite-range 
BW analogues are not necessary to be well-defined dis- 
tributions in the squeezing limit. However, in this limit, 
the corresponding total Hamiltonian can be defined as a 
self-adjoint operator. The significant difference from the 
other point interactions studied so far is the appearance 
of the cancellation of divergences in the zero-range limit, 
namely, the cancellation of singularities emerging from 
the barrier and the well terms resulting in the existence 
of discrete "allowed" values of potential strengths. 

More precisely, using the approach based on point in- 
teractions, a new type of tunneling through thin struc- 
tures, which are composed from adjacent barriers and 
wells, has been observed. As particular examples, the two 
rectangular models with the potentials aAf(x), which 
are schematically shown in Fig.[TJ have been studied both 
analytically and numerically. For each of these models, a 
pair of equations with respect to the potential strength 
constant a and the "well" parameters has been de- 
rived in the e -> limit: Eqs. JTSJ and (O for aA+(x), 
and Eqs. (ITO1) and (j2"0)) for aA~{x) permitting three dis- 
crete sets of solutions E 1 * 1 and E'. The potential aAf(x) 
is fully transparent on the set E + , whereas on the set E' 
it admits only a partial transmission as clearly demon- 
strated by Fig. [2J Contrary, due to the spatial symmetry, 



the potential aA~ (x) is fully transparent on both the sets 
E~ and E'. Outside these three sets, both the potentials 
are opaque acting as a fully reflecting wall. 

As follows from Eqs. (fT5|) - (l2T)|) . the admissible values of 
the potential strength constant a do not depend on the 
energy of an incident particle. This seems to be unusual, 
however, as illustrated by Figs. [3] and [H the resonant de- 
pendence on k appears, when the potential parameters 
are finite. For better understanding the role of adjacent 
wells in the resonant tunneling, we have introduced the 
parameters cr^ = cr, < cr < oo, into the potentials 
aAf(x), being responsible for the well depth, so that 
<7~ (cr + ) controls the well potential if a > (a < 0) as 
defined below Eqs. In the cr — > limit, the poten- 
tials aAf(x) are transformed to a typical DB structure. 
While performing this limit, the solutions to Eqs. (fTE)) - 
(f20| continuously run to infinity and Figs. [4] and [6] clearly 
illustrate the transmission difference between the cases 
a = l and a = 0. Vice versa, when one or two wells 
appear, i.e., cr starts to increase from zero, the right tails 
of the ridges shown in Fig. [5] are approaching the a-axis 
at infinity and with further increasing a, the countable 
series of these tails runs to the origin a = as a — > oo, 
so that the (a, &)-relief is becoming deformed and "ro- 
tated", taking finally the form shown in Fig. @] This 
behavior explains why it is impossible to obtain a point 
interaction from a DB potential (without wells). 

Concerning possible physical implications of the res- 
onance properties described above, one can notice that 
according to Eqs. (ITg)) -(l2Tj |) . the location of the resonance 
sets T,^ and E' is very sensitive while varying the well- 
potential part, i.e., the parameter o~ = cr if a > 
(cr + = cr if a > 0). Therefore it is reasonable in practice 
to design the barriers as system parameters given by the 
strength a, whereas the wells as adjustable external po- 
tentials controlled by the parameter cr. Inspired by the 
recent work of Turek and Cheon [17J , one can adopt the 
well part as an external attractive potential controlled by 
cr. In other words, the role of the barrier potential part 
will play a, whereas the well potential will be controlled 
by the parameter aa. Obviously, for a given a belong- 
ing to the sets E 1 * 1 or E', one can adjust a certain value 
of cr at which the barriers will be fully transparent for 
the flow of quantum particles. Any deviation from this 
value will block up the flow completely. In this way, one 
can design a quantum device like a transistor control- 
ling the flow of electrons through barriers by an external 
potential. The situation is very similar to a spectral fil- 
ter designed in Ref. [l7| where a sharp resonance termed 
the "threshold resonance" has been observed when the 
energy of incident particles equals the applied external 
potential. Roughly speaking, due to the "rotation" of 
the resonance picture illustrated by Figs. @] and [HI the 
wave number k is replaced by the strength a and this 
is a main difference between the a-resonances described 
in our case and the threshold A-resonances admitting the 
spectral filtering transmission. 

Next, similarly to the quantum Turek-Cheon device 
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(a) 




(b) 



(c) 




G 



FIG. 8: Graph with a vertex (filled ball), which describes 
a repulsive barrier junction, and three entering edges, at the 
endpoints of which terminals (empty balls) are placed where 
attractive potentials controlled by the parameters <ti , 02, and 
03 are switched on. 




(d) 



FIG. 7: Graphs that correpond to the double <5'-like point in- 



teractions obtained by squeezing the regular functions A* (a;). eters <7 
Graphs (a) and (b) describe the point interaction obtained 
from the potential aA.£(x) with a > and a < 0, respec- 
tively. Similarly, graphs (c) and (d) correspond to the poten- 
tial aAj(x) with a > and a < 0, respectively. The filled 
balls mean barriers and the empty balls together with edges 
correspond to the external potentials controlling by parame- 
ter a. 



}17| , the BW structures in which the role of the well plays 
an external potential can schematically be depicted as 
quantum graphs. In our particular case, the graphs that 
correspond to Fig. []] are shown in Fig. [7] Here, the lines 
(a), (b), and (c) describe the device with two barriers 
shown by the filled balls which are connected through 
a terminal (shown by the empty ball) with an external 
potential controlled by the parameter a. On the lines (a) 
and (b), the same potential is also applied at the side of 
one of the barriers. In fact, this is one device and the 
difference concerns its input and output as well as the 
parameter interchanging tf>i _1 . 

Of particular interest for resonant tunneling is a single 
barrier with two adjacent identical wells schematically 
shown in Fig. EJd). First, the phenomenon of full res- 
onant tunneling across a single barrier seems to be un- 
usual, but the presence of two identical wells provides a 
full transmission at certain quantized values of the exter- 
nal parameter a. Second, from a technological point of 
view such a device can be realized much easier than the 
others. A single barrier can be created simply by a nat- 
ural contact (junction) of two quantum wires, whereas 
an external potential can be applied at the endpoints 



of edges entering this junction. By the term "junction" 
we mean a part of the system from which several wires 
emerge. Varying a, one can permit the quantum flow 
to pass across this junction totally or block it up com- 
pletely. Obviously, this system with two terminals may 
be generalized to several terminals like that schemati- 
cally illustrated by Fig. |8j where in general the param- 



3 > 



1, 2, 3, are different and any of them may 
be assumed zero (the absence of wells). 

Thus, we have treated the two simplest examples of 
multilayer structures which are illustrated by Fig. Q] and 
in the zero-range limit can be approximated by point in- 
teractions with total resonant tunneling. The tunneling 
effect observed in this paper by using the point interac- 
tion approach may be used in more realistic structures 
like a resonant tunneling diode with an additional well 
termed "prewell" [3l[ . Graph (b) in Fig. [7] corresponds 
to this device if the external potential is applied only 
to the edge in front of the DB structure. Nevertheless, 
the structures with the well located between the barriers 
should be of interest as well because in this case one can 
get a full resonant tunneling. In general, the external 
potential parameters <Jj, j = 1, 2, ... , can be assumed 
different and similar methods can be used for studying 
more complicated geometries. We hope to pursue this 
point in the future. 

The phenomenon described in this paper should be 
observable not only in quantum systems, but also in sys- 
tems in which classical waves can propagate. In view 
of the formal analogy between the Schrodinger equation 
and the electromagnetic Helmholtz equation, this phe- 
nomenon can be directly applied to a medium with one- 
dimensional stepwise inhomogeneity of dielectric perme- 
ability. While the experiments involving electrons are 
usually difficult to realize (mainly due to the smallness 
of the electron de Broglie wavelength at usual temper- 
atures), the observations on microwave tunneling pro- 
vide clear data. Besides this, for instance, for TE waves 
(E x = H y = E z = 0, an incident wave and inhomogene- 



9 



ity coincide with the x direction), the analogue of the 
strength a appears to be the frequency of electromag- 
netic field which can be tuned (32j |. 
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